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ABSTRACT 

We have computed models for ultraluminous X-ray sources ("ULXs") consisting of a black-hole accretor of 
intermediate mass ("IMBH"; e.g., ~1000 Mq) and a captured donor star. For each of four different sets of 
initial donor masses and orbital separations, we computed 30,000 binary evolution models using a full Henyey 
stellar evolution code. To our knowledge this is the first time that a population of X-ray binaries this large has 
been carried out with other than approximation methods, and it serves to demonstrate the feasibility of this 
approach to large-scale population studies of mass-transfer binaries. In the present study, we find that in order 
to have a plausible efficiency for producing active ULX systems with IMBHs having luminosities > lO"*'^ ergs 
s"', there are two basic requirements for the capture of companion/donor stars. First, the donor stars should be 
massive, i.e., > 8 Mq. Second, the initial orbital separations, after circularization, should be close, i.e., < 6-30 
times the radius of the donor star when on the main sequence. Even under these optimistic conditions, we show 
that the production rate of IMBH-ULX systems may fall short of the observed values by factors of 10-100. 

Subject headings: (stars:) binaries; general — X-rays: binaries — galaxies: star clusters — accretion, accretion 
disks — black hole physics 
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1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are off-nucleus 
sources in external galaxies with luminosities of 10^^ < < 
10^' ergs s~'. They have been discovered in great numbers 
with ROSAT, Chandra, and XMM-Newton (Fabbiano 1989; 
Roberts & Warwick 2000; Ptak & Colbert 2004; Fabbiano & 
White 2004; Colbert & Miller 2004). For reference, the Ed- 
dington limit of a black hole accretor of ~10 Mq is in the 
range of ~ 1.5-3 x 10'''' ergs s"', depending on the chemical 
composition of the accreted material. ULXs are especially 
prevalent in galaxies with starburst activity, including ones 
that have likely undergone a recent dynamical encounter (e.g., 
Fabbiano, Zezas, & Murray 2001; Wolter & Trinchieri 2004; 
Fabbiano & White 2004; Colbert & Miller 2004). Many, but 
not all, of the ULXs are identified with young star-forming 
regions (see, e.g., Fabbiano & White 2004). 

The central question regarding this important class of 
sources is whether they represent an extension in the lumi- 
nosity function of binary X-ray sources containing neutron 
stars and stellar-mass black holes (BHs), or a new class of ob- 
jects, e.g., systems containing intermediate-mass black holes 
(10^- 10* Mq- hereafter "IMBHs"). 

Colbert & Mushotzky (1999) first suggested that some 
ULXs harbor IMBHs. The motivation for this is clear. The 
Eddington limit for a 1000 Mq black hole is ^ 10*^ ergs s"', 
compared to only ^lO^*** ergs s"' for a typical neutron star, 
and '^lO times this value for a stellar-mass black hole accre- 
tor. Moreover, the spectra from IMBHs might be expected to 
have low inner-disk temperatures, as is inferred for some of 
the ULXs (Miller et al. 2003; Miller et al. 2004; Cropper et 
al. 2004). 

A number of ideas have been put forth for ways to circum- 
vent the problem of how ^10 Mq BHs could have apparent 
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L, values as high as lO'*^- 10'" ergs s"'. King et al. (2001) 
suggested that the radiation may be geometrically beamed 
by a thick accretion disk so that the true value of does 
not, in fact, exceed the Eddington limit. Kording, Falcke, & 
Markoff (2002) proposed that the apparently super-Eddington 
ULXs are actually emission from microblazar jets that are rel- 
ativistically beamed along our line of sight. However, studies 
of the giant ionization nebulae surrounding a number of the 
ULXs (Pakull & Mirioni 2003) seem to confirm the full lu- 
minosity inferred from the X-ray measurements. Begelman 
(2002) and Ruszkowski & Begelman (2003) found that in ra- 
diation pressure dominated accretion disks super-Eddington 
accretion rates of a factor of ~10 can be achieved due to the 
existence of a photon-bubble instability in magnetically con- 
strained plasmas. 

In two previous studies we have explored the evolution of 
binary stars consisting of 2- 17 Mq donor stars and stellar- 
mass BHs as possible model systems for ULX sources (Pod- 
siadlowski, Rappaport, & Han 2003; Rappaport, Podsiad- 
lowski, & Pfahl 2005a). These BH binaries are assumed to 
form from massive primordial binaries in regions of intense 
star formation. In our earlier work we showed that the forma- 
tion efficiencies for such systems were adequate to explain the 
observed formation rates of ULX systems. However, the Ed- 
dington limit would somehow have to be violated by a factor 
of ~10-30, and even then, the mass transfer rates were only 
marginally sufficient to explain the luminosities of the most 
powerful ULXs. 

In this work we investigate the evolution of a population 
of binaries containing an IMBH with a representative mass 
of 1OOOM0. The mechanisms by which an IMBH is created 
in star-forming regions and acquires a companion/donor star 
are highly uncertain. One proposed scenario for their forma- 
tion in star clusters involves a runaway stellar collision pro- 
cess leading to the production of a very massive star which 
somehow evolves to core collapse before much of its envelope 
has been lost. While the dynamics of the runaway collision 
may be understood quantitatively (e.g., Portegies Zwart & 
McMillan 2002; Portegies Zwart, Dewi, & Maccarone 2004; 
Giirkan, Freitag, & Rasio 2004) there is a great deal of uncer- 
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tainty over how, and whether, such a massive star evolves to 
core collapse and the formation of an IMBH. If we allow for 
the presence of an IMBH at the center of a star cluster shortly 
after the formation of the cluster itself, then one must con- 
sider how the IMBH acquires a mass-transferring companion 
star. Presumably the IMBH could capture a companion from 
stars near the cluster core via either tidal capture or exchange 
encounters with primordial binaries in the cluster (Hopman, 
Portegies Zwart, & Alexander 2004; Pfahl 2005). After the 
IMBH has acquired a companion, but before steady mass 
transfer is possible, the star must survive tidal circulariza- 
tion and, more specifically, the tidal heating that accompanies 
the process. There has been much discussion in the literature 
about this issue (e.g., Ray, Kembhavi, & Antia 1987; Podsiad- 
lowski 1996; Alexander & Morris 2003). Finally, the binary 
system has to survive encounters with other passing stars — 
either single or binary stars. For earlier work on IMBH-ULX 
binary evolution see, e.g., Rappaport, Podsiadlowski, & Pfahl 
(2005b), Kalogera et al. (2005), and Patruno et al. (2005). 
For a very recent investigation of the dynamical capture of 
companion/donor stars by an IMBH, see Blecha et al. (2005). 

Blecha et al. (2005) found a wide range of semimajor axes 
for the companion stars captured by the IMBH, with the ma- 
jority in the range ^ 100- W^Rq (and ^90% lying between 
30 and 3 x 10"* Rq). We know from our previous work, and 
will show in the present study, that very wide binary systems 
do not produce ULXs with interesting durations during the 
lifetime of the star cluster Therefore, we have chosen to 
study systems with initial orbital separations in the range of 
- 20- 1000 wherein ~40% of the Blecha et al. (2005) 
systems are found. We further assume that the IMBH has 
only a single stellar companion which is counter to the main 
findings of Blecha et al. (2005). However, the assumption of 
a single companion would be quite reasonable for Model C 
(low stellar density) of Blecha et al. (2005) and for times up 
to ~30 Myr, especially when taking into account that these 
authors assumed a primordial binary fraction of 100%. 

In this work we follow 30,000 binary evolution sequences 
for each of four different distributions of starting masses and 
orbital periods. The close orbits are assumed to have circu- 
larized by tidal friction. Each binary is followed with a full 
Henyey stellar evolution code until the entire envelope of the 
donor star has been transferred. To the best of our knowledge, 
this is the first time that a binary population study this large 
has been carried out with a full Henyey code. ^ 

2. POPULATION STUDY 

The stellar evolution of the donor stars, including mass 
loss, was followed with E Z which is a stripped down, rewrit- 
ten version of a subset of the stellar evolution code devel- 
oped by P. P. Eggleton (Paxton 2004). The physics of the 
program is unchanged from Eggleton's (essentially as de- 
scribed in Pols et al. 1995), but the structure of the code 
has been modified to facilitate experiments involving pro- 
grammed control of parameters. There are zero-age main- 
sequence (ZAMS) starting models for a variety of metallici- 
ties (from Z = lO""* to Z = 0.03) and masses (from 0.1 to 100 
Mq), with arbitrary starting masses created by interpolation. 
A user-provided procedure is called between steps of the evo- 
lution to inspect the current state, to make changes in param- 
eters, and to decide when and what to record to log files. The 

^ Prior to this, the largest number of mass-transfer binaries computed in a 
single study (in the published record) was 5500 (Nelson & Eggleton 2001). 



TABLE 1 
ULX-IMBH Population Models 



Model 


A/low'' 




/max 




A 


5 


50 


15 


40-900 


B 


1 


10 


15 


20-600 


C 


5 


50 


3 


40-200 


D 


1 


10 


3 


20-130 



"Lower limit on the donor mass, in units of Mq 
''Upper hmit on the donor mass, in units of Mq 

■^Maximun initial orbital sepai'ation in units of that required for the donor 
star to fill its Roche lobe; / = fli„it/oRL 
''Range of initial orbital separations, in units of Rq 



source code and data for E Z can be downloaded from the web 
at <http://theory.kitp. ucsb.edu/^paxton>. For all models in 
this particular study, the number of stellar mesh-points was 
fixed at 200 in the interest of minimizing computation time. 

A binary "driver" which computes the changes in orbital 
separation and the mass loss and transfer rate (according to 
the appropriate conservations laws) was written specifically 
for this project. However, it is similar to the procedures and 
approaches which we have used in our stellar evolution pro- 
grams in the past (see, e.g., Podsiadlowski, Rappaport, & 
Pfahl 2002). In particular, we utilize a semi-explicit scheme 
for determining the mass loss at each time step based on the 
instantaneous location of the Roche lobe within the atmo- 
sphere of the donor star, which is approximated as having an 
exponential density distribution. 

For each binary we choose, via Monte Carlo means, the ini- 
tial donor mass uniformly between Mio„ and Mhigh (see Table 
1 for specific model values). The initial orbital separation was 
chosen uniformly between 1 and /^ax times the separation re- 
quired for a ZAMS star (of the same mass) to fill its Roche 
lobe (see Table 1). For Models B and D the donor masses 
ranged between 1 and IQMq, while Models A and C featured 
much higher mass stars (5-5OM0). Models A and B had 
initial orbital separations up to /nj^x = 15, which crudely rep- 
resent what might be expected from the tighter exchange en- 
counters with passing primordial binaries (see also Blecha et 
al. 2005). In contrast. Models C and D had /max = 3 to roughly 
approximate the case of direct tidal capture, and circulariza- 
tion, of a single field star by the IMBH (see, e.g., Hopman 
et al. 2004). Parameters for the four illustrative distributions 
we have chosen to explore are summarized in Table 1. We 
ran 30,000 binary evolution sequences for each of the models 
listed in Table 1 . 

All models were run using 60 nodes of the elix3 Beowulf 
cluster located at the University of Sherbrooke, Quebec. The 
run time for each of the 4 models was ~30 hours. 

3. IMBH-ULX RESULTS 

For each evolution step of each IMBH-ULX binary, the X- 
ray luminosity is computed according to L^■ = rjMc^, where 
M is the mass transfer rate, and rj is the BH conversion ef- 
ficiency from rest mass to radiant energy, computed accord- 
ing to the spin of the black-hole. All IMBHs are assumed to 
start with zero angular momentum. We then follow the spin 
history of the IMBH (according to the formalism of Bardeen 
1970); however, because of the extreme mass ratio between 
the IMBH and the donor star, the dimensionless spin parame- 
ter of the IMBH cannot exceed ^0.07, and the corresponding 
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Log Evolution Time (yr) Log Evolution Time (yr) 

Fig. 1 . — Simulated evolution of the X-ray luminosity function of IMBH black hole binaries with time since the birth of the host star cluster. For each panel, 
evolution tracks from 30, 000 X-ray binaries were computed and then registered in each of the 700 X 700 pixels that are traversed. In the calculations used to 
produce these plots, the Eddington limit (at 2 X 10'*' ergs s"') was enforced; however, for clarity the L, values displayed here are without any Eddington limit. 
This affects only a relatively small fraction of the highest luminosity systems. The colors represent the square root of the relative populations, with red through 
red corresponding to actual ratios of ~100 to 1. The distributions of initial masses and orbital separations for the four models are specified in Table 1. 



value of T] does not grow appreciably above that expected for 
a Schwarzschild BH (1 - \/8/3). 

For each evolution step of each IMBH-ULX binary, the X- 
ray luminosity and corresponding evolution time are recorded 
in a 700 x 700 matrix, by simply adding "1" when the track 
crosses a particular array element. The contributions from all 
of the tracks are added together. The matrix covers 7 decades 
in and 3.5 decades in time, in equally spaced logarithmic 
intervals. The results are shown in Fig. 1 as color images 
with a square-root scaling in intensity to enhance the dynamic 
range. Time zero is arbitrarily taken to be a common forma- 
tion epoch for all stars in the cluster, as well as the capture 
time of a companion by the IMBH. 

Figure 2 shows more quantitative plots of the results for 
the models "imaged" in Fig. 1. The curves in each plot show 
the numbers of IMBH systems that are transferring mass at 
any given time - after their formation - with values of L.i that 
exceed lO-'^, lO-'^, and 10'*" ergs s~'. The initial orbital peri- 

3/2 

ods of the systems being computed are typically ~ 1 /max 
days, and end with Porb often exceeding a year Note that even 
for Models A and C, which emphasize massive donor stars, a 
maximum of ^4% and ^^25%, respectively, of systems have 
Lx > 10'*° ergs s"', but only for a brief interval of ~1 Myr 



during the star cluster's active lifetime. 

From Figs. 1 and 2 we see that Models A and C, with 
higher-mass donor stars, produce generally higher values 
than Models B and D, and the systems are most "active" be- 
tween 3 and 30 Myr By contrast. Models B and D, with 
lower-mass donors, typically yield significantly lower values 
of Lv, but these systems appear over a much longer range of 
timescales, i.e., out to ^ lO'' yr. For the latter models, there 
are ULX-luminosity sources between ^50-500 Myr due to 
donor stars that are ascending the giant branch; however, these 
are probably too rare (fewer than ~0.003 of all the systems) 
to explain many ULXs. These differences between systems 
with higher and lower mass donor stars can be explained by 
the fact that higher mass stars have shorter nuclear evolu- 
tion timescales. This leads to shorter overall lifetimes of the 
donors and correspondingly higher values of M. Systems on 
the giant branch, in all four panels, but especially Models B 
and D, can be recognized as a red band above the main body 
of systems in which the donor star is still on the main se- 
quence. The vertical edge at ^^23 Myr in Model B (Fig. 1) 
marks the evolution time required for a 10 Mq star (the most 
massive in the model) to ascend the giant branch. Due to the 
wide initial separations, most of the systems in Model B do 
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Log Evolution Time (yr) 

Fig. 2. — The numbers of IMBH binaries as a function of time for three 
different lower limits on the X-ray luminosity. Blue and purple curves - all 
systems. Green and cyan curves - Lx > 10^' ergs sec"' . Red and orange 
curves - L, > 10'"' ergs sec"' Values are derived from the results in Fig. 1. 
Color scheme: Models A and C (blue-green-red); Models B and D (purple- 
cyan-orange). In each case, the number of active systems at any given time, 
is divided by 30,000 (the number of binaries evolved) to provide an estimate 
of the fractional "efficiency" for producing IMBH-ULX binaries. 



not commence mass transfer until the donor evolves well be- 
yond the main sequence. 

The distribution of initial orbital separations is equally im- 
portant in determining the production of ULXs as is the donor 
mass. Models C and D, with closer initial orbital separations, 
have generally significantly larger populations of active ULXs 
than for their counterpart Models A and B. This is true since in 
initially wide systems - by the time mass transfer commences 
- the donor star is already quite evolved, and the remaining 
evolution time on the giant branch is rather short. Hence, 
these systems do not persist as ULXs for very long, and there 
is a correspondingly small probability of finding the system 
in an active state. 

Finally, we comment on a technical issue of the binary evo- 
lutions. We find that for some binary models the donor star 
envelope is not completely transferred to the IMBH before 
the stellar code terminates. This occurs predominantly for the 
most massive donor stars (i.e., initial masses > 4-QMq) and 
wide initial orbital separations (i.e., / > 7, where / is defined 
in Table 1); thus, it affects mostly Model A. The computa- 
tional problems begin to occur close to the point of core He ig- 
nition when the rapid readjustment of the internal stellar struc- 
ture, coupled with high mass transfer rates, require very small 
time steps. This problem is not insurmountable, but one that 
may require significant additional computational time. Over- 
all, we estimate that we have successfully accounted for the 
following fractions of ULX-luminosity systems in Models A, 
B, C, and D: 66%, 96%, 92%, and 90%, respectively. Thus, 
the appropriate corrections could be made for the correspond- 
ing ULX production rates that we present. 

4. SUMMARY AND CONCLUSIONS 

In this work we have explored a number of different models 
for ULX sources with intermediate-mass BH accretors. These 



involved four different ranges of donor masses and initial or- 
bital periods. For each model, the evolution of a very large 
number of individual binary systems was computed. We find 
that in order to have a plausible efficiency for producing active 
ULX systems, there are two basic requirements for the cap- 
ture of companion/donor stars. First, the donor stars should 
be massive, i.e., > 8 Mq. Second, the initial orbital separa- 
tions, after circularization, should be close, i.e., ~6-30 times 
the radius of the donor star when on the ZAMS. If these highly 
favorable conditions do occur, we see from Fig. 2 that up to 
-25% of the IMBHs could have > 10"*" ergs s"' , but only 
for a brief interval of time on the order of 1 Myr. This is 
at least qualitatively consistent with the fact that such objects 
are quite rare, occurring, on average, only once in every —100 
galaxies (see, e.g., Ptak & Colbert 2004). 

We make the following crude estimate of the number of 
ULXs with L^ > 10'^^ ergs s"' that might be found in a typical 
spiral galaxy: 

MmBH-ULX — A'ysc • /iMBH ' /cap ' /on ' /dur , ( 1 ) 

where A^ysc is the steady-state number of young (i.e., <50 
Myr), massive star clusters potentially capable of producing 
an IMBH, /mbh is the fraction of all such clusters which pro- 
duce an intermediate-mass BH, /cap is the fraction of IMBHs 
that capture a massive star into a close orbit, /□„ is the frac- 
tion of systems that would be "on" with this luminosity at a 
particular time (see Fig. 2), and /dur is the fraction of the 
cluster lifetime when the massive stars would be transferring 
mass to the IMBH at a sufficiently high rate (see Fig. 2). 
For illustrative parameter values we might take A'ysc — 10 (if 
anything, a generous estimate), /mbh — 0.1, /cap — 0.05 (see 
Blecha et al. 2005), and /□„ x /dur ranges between 0.001 
and 0.02 (as derived from Fig. 2). We would then find 
Mmbh-ulx — 5 X 10"^- lO""'. Thus, only one normal galaxy in 

- 10-^- 10^ might be expected to harbor a ULX with L, > 10"*" 
ergs s"', about a factor of 10-100 lower than is indicated by 
the observations (see, e.g., Ptak & Colbert 2004). However, 
needless to say, a number of these parameter values are highly 
uncertain. 

We have not considered in any depth either the formation of 
the IMBH itself, or the capture of the companion/donor star. 
Nor have we taken into account the possibility that more than 
one companion star could be captured by the IMBH within 
the same time interval, and these could interact dynamically 
(see Pfahl 2005; Blecha et al. 2005). However, at least for (i) 
modestly low cluster central star densities (— 3 x 10"* pc"^*), 
(ii) the duration of the evolutionary timescales of the massive 
donor stars of interest, and (iii) a reasonable binary fraction 
of —50%, this may be a good approximation. In summary, 
we conclude that in order for the IMBH model of ULXs to be 
viable, at a minimum, there must be a high efficiency for pro- 
ducing IMBHs in young star clusters, and a high probability 
of capturing a massive companion star into a relatively close 
orbit. Even then, however, it is not obvious that the IMBH- 
ULX model can account for the observed systems. 

Finally, even for our models with massive donor stars, 
IMBH systems with 10^^ < L, < lO'*" ergs s"' are only 

— 1.5-2.5 times more numerous than systems withLj > 10"*" 
ergs s"'. This ratio is inconsistent with the rapid falloff in 
ULXs with luminosity that is observed. 

We thank Alan Levine for helpful discussions, and Robert 
Harris for his participation at an early stage of this study. 
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